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Abstract 

We present an analysis of the quantum fluctuations around the electroweak 
sphaleron and calculate the associated determinant which gives the 1— loop 
correction to the sphaleron transition rate. The calculation difl"ers in various 
technical aspects from a previous analysis by Carson et al. so that it can be 

considered as independent. The numerical results difl"er also — by several 
orders of magnitude — from those of this previous analysis; we find that the 
sphaleron transition rate is much less suppressed than found previously. 
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1 Introduction 



The electroweak theory is known |^ for quite some time to possess a topologically 
nontrivial solution which describes a saddlepoint between two topologically distinct vacua. 
The recent interest in this solution has been centered around its possible role in generating 
baryon number violating processes in the early universe or even at accelerator energies 
0-0]. The rate of sphaleron transitions in the range of temperatures Mw(T) <^ T <^ 
MwiT)/au] has been derived on the basis of the work of Langer P and Affleck 0] by 
Arnold and McLerran It is given by 

r = ^J^e-^'^'^^K. (1.1) 
2n 

Here U- is the absolute value of the eigenvalue of the unstable mode, the prefactor Af refers 
to normalizations introduced by the translation and rotation zero modes and is given in 
detail below. Ed is the classical sphaleron energy and the factor k, takes account of the 
quantum fluctuations of the sphaleron. It is given by 

detAfpdetA^ /2 

MetA^detA^^y ^ ^ 

where the symbols A denote the small fluctuation operators. They are obtained by ex- 
panding the gauge fixed action (gf) and the Fadeev-Popov action (FP) evaluated around 
the sphaleron (S) and the vacuum (0), respectively. If the fluctuation operators are diago- 
nalized the determinants are formally given by the product of the squared eigenfrequencies 
A^. The determinant det A^j of the gauge fixed action is to be evaluated without the zero 
modes and with the eigenvalue of the unstable mode, = — tu^ replaced by its absolute 
value. 

The first evaluation of k, by Akiba, Kikuchi and Yanagida was restricted to the 
three lowest partial waves. The authors concluded from their results that the inclusion of 
this correction did not suppress the sphaleron transition rate. Of course they considered 
their conclusion as tentative and to be checked by a more precise evaluation. Subsequently 
Carson and McLerran evaluated k in an approximation scheme (referred to as DPY in 
the following) developed by Diakonov, Petrov and Yung |T0[, a complete exact numerical 
evaluation was presented first by Carson et al. [|lT| (referred to as CLMW in the following). 
The results of this exact calculation differ significantly from the DPY approximation and 
from a perturbative estimate. It is therefore of interest to repeat this exact evaluation 
and this is - besides a general analysis of the fluctuation Lagrangean and its partial wave 
reduction - the subject of our investigation. 

Our investigation departs from the one of CLMW in various points. While we use the 
background gauge as CLMW, we use another angular momentum basis and a different 
scheme for evaluating the determinant. 

The analysis of the small fluctuations around the sphaleron requires a partial wave 
decomposition with respect to the quantum number K = J + I (i^-spin). Our present 
work is based on the analysis of Ref. |TB| where the sphaleron stability was investigated. 



1 



The small fluctuation equations obtained there had to be modified however. While that 
investigation was performed in the Aq = gauge we use here the background gauge. So 
a gauge fixing term had to be added and also we had to construct the small fiuctuation 
Lagrangean for the Fadeev-Popov modes. The evaluation of the determinant has been 
performed using the Euclidean Green function technique in analogy to some recent inves- 



tigations of one of the authors (J.B.) |T^. The results of our calculation have already been 



communicated previously in a short version |T^. Here we present an extensive version 
including the complete partial wave analysis. 

The paper is organized as follows. We present the basic relations for the electroweak 
theory and the sphaleron solution in section 2. The small fiuctuation expansion and its 
partial wave analysis are discussed in sections 3 and 4, respectively. The explicit equations 
of motion which constitute our first main result are presented in Appendix A. In section 
5 we collect the main formulae for the 1-loop correction to the sphaleron transition rate 
and define more precisely the zero mode prefactors and the fiuctuation determinant, some 
material being deferred to Appendix B. In section 6 we express the fiuctuation determinant 
in terms of Euclidean Green functions, a formulation that is used as the basis of our 
numerical evaluation. Sections 7 and 8 contain the discussion of renormalization and of 
the treatment of the zero and unstable modes. Some aspects of the numerical calculations 
and the results are presented in section 9 and in Appendix C. 



2 Basic Relations 

The action of a pure SU(2) gauge theory with minimal Higgs sector is given in Minkowski 
space as 

S = J d^xC{x), (2.1) 
where the Lagrangean density is given by 

We have used a bar to denote the original fields and coordinates (x, $, . . .), we wfll 
use the same letters without bar (x, F^^, . . .) for a rescaled version of these quantities 
(see below). 

The non-Abelian field strenght tensor and the covariant derivative of the Higgs field 
are given by 

and V, = d^-'-gr^Wl (2.3) 

respectively, where are the Pauli matrices. 

After spontaneous symmetry breaking the vacuum expectation value of $ is f 7^2 and 
the W^^ bosons get a mass = gvj^. The mass of the physical Higgs boson becomes 
mil = V2Xv. We define further, as usual, the weak fine structure constant as = g'^/4:7r. 
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In view of the renormahzation of the theory at finite temperature it is suitable to rescale 
fields and coordinates as 

As a result of this rescaling the action takes the form 

S= \ f d*xC{x) (2.5) 
g^ J 

with the Lagrangean density 

C = -^FI^^F^"^ + 2 (P^$)^ (P'^*) - ($^$ - l)^ , (2.6) 

Here we have introduced the ratio of Higgs and W masses ^ via 

mjj _ 8A 
rriw g 

Field strength and covariant derivative reduce to 

-D, = (2.8) 

As a consequence of this rescaling the action does not depend any more on the vacuum 
expectation value of the Higgs field and the couplings appear only in the mass ratio ^. As 
explained in p, |lT| at high temperatures the fields can be considered as essentially static 
and the theory reduces to a three dimensional theory if the vacuum expectation value is 
rescaled to the temperature dependent one [llB, M 



v{T) = v{0)Jl-TyT,\ (2.9) 



where = 2v{0f 



1 + ^ 
8A 



(2.10) 



is the critical temperature for symmetry restoration. While the masses become then func- 
tions of the temperature, Mw{T) and Mh{T), this dependence cancels in their ratio ^. 
One obtains for the three dimensional Euclidean action 



Se = - I d'xCEix), (2.11) 

gs 



with the Lagrangean density 



^Eix) = \F;,F,,a + 2 (D^$)t (P^$) + ($^<^> - l)' • (2.12) 



^Our scale units differ from those of CLMW, we use mw instead of gv as tfie basic unit. 
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The effective coupling constant of the three dimensional theory which describes the 
interaction at high temperatures is given by 

For temperatures 

T < Mw{T)/a^ (2.14) 

this coupling is much smaller than 1 and an expansion with respect to should be reliable. 

The variation of the action ( pl5| ) with respect to VF^ and leads to the classical 
equations of motion 







and V^'V^^ + |^$t$ _ $ = q, (2.15) 

where the covariant divergence of the field strength tensor is given by 

(P-F^^)" = S'^f;, + e^^W^'F;, (2. 16) 

In a static saddle point solution to these equations has been constructed explicitly, 
the well known sphaleron. We choose here an Ansatz that differs from the one of Ref. 
by a SU{2) rotation of the Higgs field by an angle 7r/2 (which is another special case of a 
more general parametrisation |^) : 

= 

fA{r)-l 

$ = H,{r){^\ (2.17) 



1, 



with 



Xm = — . (2.18) 
r 

Our solution will then be a gauge rotated version of the one given in 0. 

Inserting this Ansatz into the classical equations of motion ( |2.15| ) one obtains for the 
profile functions /a and Hq a coupled system of the form 



/a(/1 - 1) 



r 
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< = leHo{H'o-l)-lH',+ ^-l^^Ho. (2.19) 



The energy of the sphaleron configuration is given by 



E,, = ^^ IdrHir), (2.20) 



whith the Hamiltonian density 



n{r) = r\ + ^(/l - 1)^ + 2r'H'l + {U - l?Hl + \er\Hl - if. (2.21) 



Requiring the finiteness of the energy yields the following boundary conditions for the 
profile functions 

/a(0) = -1 
/a(oo) = 1 
i^o(O) = 

Hq{oo) = 1. (2.22) 

The numerical solutions of the equations ( p.l9| ) used later on were obtained using a 
method developed in |]18| . The classical sphaleron energy as a function of C, is presented in 
Fig. 1 in units of Mw{T)/aw The classical Euclidean high temperature action as defined 
in Eq. ( |2.11|) is then given by 

5| = E,i/T. (2.23) 



3 Small Fluctuations around the Sphaleron 

In this section we will develop the small fluctuation expansion around the sphaleron. We 
will start with expanding the Lagrangian up to second order in the small fluctuations and 
fix the gauge. Subsequently we will expand the small fluctuations into partial waves. 
The gauge and Higgs fields of the theory are expanded around the spaleron configuration 

via 

W^o = < 

Wt = At + at (3.1) 

= + (3.2) 

und $c« denote here the classical solution ( |2.17| ) 

Aa _ fA - 1 . 
r 

$d = Ho(^^. (3.3) 

and U{(f)) is given by 

U{ip) = exp(irV^). (3.4) 
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So the Higgs field fluctuations are parametrized by the isosinglet h and the isotriplet (p^ 
and those of the gauge field by a^. We note that $ has to be expanded up to second order 
in the fluctuation fields if the second order Lagrangean is to be determined; so includes 
second order terms ^. 

Inserting these expansions into Eq. (|2.6| ) and collecting terms of the same order in 
the small fluctuations we find in zeroth order the classical sphaleron action, the first order 
contribution vanishes since the sphaleron configuration is a saddle point of the action. The 
second order Lagrangean density in which we are interested here becomes 



(3.5) 



The tildes on the covariant derivatives indicate that only the classical gauge fields are to 
be used: 



pa 



d^Al-d^A^ + e'^'^AlAl 
d,al + e^'^Alal. 



(3.6) 



In order to eliminate the gauge degrees of freedom we have used the background gauge 
which restricts the fluctuating fields. In general form the three constraints are given by 







(3.7) 



for a = 1,2,3 where ^^^^ is defined in Eq ( p.2| ). With our parametrization of the Higgs 
field they take the form 

r''={v^a^'Y + 2Hop\ (3.8) 

Here only the linear terms in ^'^^^ had to be taken into account since only the square of JF^ 
appears in the gauge- fixed Lagrangean (see below). 

The Fadeev-Popov determinant for this gauge condition is given by 



detAf,^ 



det 
det 



ab 



+ Hi 



(3.9) 



Finally one obtains the gauge fixed Lagrangean density 



(2) 

gf 



cm 



■^This has been done in exactly the same way injT^, but wrongly stated there in Eqs. (3.1) and (3.2) 
which should be replaced by our Eqs. (^^) and (4.2). 
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2 



a 



^ TT2„a„a 



1 

+AhA1diip'' - 2HQhA1al 



(3.10) 



The background gauge has led to a decouphng of the time components of the gauge 
fields. Their contribution to the Lagrangean density 



r(2) 



2 ° 



ab 



(3.11) 



leads to a fluctuation operator identical to that of the Fadeev-Popov ghost fields. 



4 Partial Wave Expansion of the Fluctuations 

In order to arrive at a suitable basis for our numerical computations we have to decompose 
our system of small fluctuations into partial waves with respect to the K spin K = J + 1 
combining angular momentum and isospin, which is conserved on the sphaleron back- 
ground. We avoid a large amount of Clebsch-Gordan algebra by chosing a basis of tensor 
spherical harmonics which is constructed using cartesian vector operators (see e.g. 
We define the following dimensionless operators 



t^. =r\7 -A =A-rV =AA 

These can be used to obtain a suitable basis of tensor spherical harmonics with the spin- 
isospin transformation properties appropriate to the different fields: 

3 

«0a = H H Sa(r)J"yft-M(5') 

_fs:,Ma=l 
9 

«ia = ^ta(r)IJ„YKAf(£) 
K,M a=l 
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KM 



= EEP"WJa>^/rM(x)L . (4.2) 
K,M a=l \^/ 



The algebraic properties of the operators J and I and of the tensor basis are discussed 



extensively in [20| and will not be repeated here. All manipulations have been performed 



using REDUCE, starting with the basic Lagrangean density ( Um) and using the substitu- 



tions ( |3.1| ), ( |3.2| ) and ( |4.2|) and the gauge fixing based on (|37 

The Lagrangean decomposes into contributions for each K spin. Varying these La- 
grangeans one obtains a set of 16 linear coupled differential equations for the radial func- 
tions Sa, ta, hi and pa- According to their parity transformation properties the 16 am- 
plitudes fall into two groups which do not couple between each other, an electric and a 
magnetic sector. The electric sector consists of the amplitudes t^, t^, tj, tg, p^, hi and S3 
with parity tt = (—1)^+^ and the magnetic sector contains the amplitudes ti, t2, t^, tg, ^9, 
Pi, P2, si and S2 with parity vr = (— 1)"^. Furthermore the ghost amplitudes si, S2 and S3 
are not coupled to the other fields. They form the Fadeev-Popov sector. So altogether we 
have a 6 X 6 and a 7 x 7 subsystem for the gauge and Higgs fields. Only si and S2 form a 
coupled system in the Fadeev-Popov sector but we will not divide this sector furthermore. 
We will refer to the three sectors by a symbol a which can have the values E, M and FP 
and the different coupled channnels will be referred to by the two quantities a for the fields 
and K for the partial wave. 

The K spins K = and K = 1 have to be considered separately: In the case K = 
the J- and I- operators act on constant spherical harmonic Iqo so that only the operators 
I^, I^, and contribute; the amplitudes t2, t^, tg, ^7, ts, tg, p2, ps, S2 and S3 have to be 
discarded and we have only 6 fluctuation equations. In the case K = 1 the action of the 
operators and on the spherical harmonics Yim oc x yields identical tensors and the 
same is true for and I^: 



^ia^j ^ia^j ^ia^j ~t~ ^ij^a ^X^X^X^ 

One of the amplitudes ^3 and and one of the amplitudes und tg have to be discarded. 
We have chosen 

tg = tg = 0, for K=l. (4.4) 

This elimination cannot be done simply in the general fluctuation equation but has to be 
performed already in the Lagrangean. 

For the following general developments we will design by a column vector which is 
formed by the amphtudes of one of the coupled sectors. The fluctuation equations can 
then be written in the general form 



^. + + _ m^^. = v,^q,._ (4.5) 



r 
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They are given explicitly in Appendix ^ As one can see there the fluctuation equations 
have some shortcomings: the amplitudes do not display a unique centrifugal term, the 
potential Vij is not symmetric and in the vacuum, i.e. for (/a 1, Hq ^ 1, ^ 0, Hq ^ 
0) the amplitudes are still coupled. One has to find therefore a transformation 

^^ = C,,^j (4.6) 

of the fields that brings them to a form 

r 

where the potential satisfies 

Such a transform must obviously exist since the Lagrangean (|3.10|) is a symmetric bilinear 
form in the fields. 

The manipulations leading to this transform have been performed using REDUCE. The 
new amplitudes, the associated angular momenta and masses and the potential in the new 
basis are also given in Appendix The equations in this form constitute the basis of our 
further developments. 



(4.7) 



und lim Vi^ = 0. 

r— >oo ■' 



(4.8) 



5 The Fluctuation Determinant k 

The rate for sphaleron transitions - based on the general theory of Langer - has been 
obtained by various authors ([0, |[ |]) to be given by 

r/V = ^T-''MtrKotVroteM-Eci/T)K. (5.1) 

ZTT 

Here n and the other prefactors arise from the Gaussian functional integration of the 
quadratic fiuctuations around the classical saddle point solution, n is essentially (see 
below) the determinant corresponding to the fiuctuation operator 



(5.2) 



The calculation of k is, besides the general partial wave analysis, the main subject of this 
publication. 

The prefactors labelled tr and rot are related to the existence of six zero modes due to 
the invariance of the theory with respect to translations and rotations which is broken by 
the sphaleron solution. These modes satisfy 

A,,^Y' = = 0- (5-3) 
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The functional integration of these modes has to be treated separately (see |, ^ p^l). 
They have to be excluded in evaluating the fluctuation determinant. These pref actors and 
their calculation are discussed in Appendix B. The results are displayed in Fig. 2. They 
agree within a few percent with those of CLMW apart from factors 2^^/^ arising from the 
different scale factors Mw resp. gv used to make the radial variable dimensionless (see Eq. 

As obvious from the discussion in the previous section K spin and parity invariance lead 
to a decomposition of the fluctuation operator into a direct sum of fluctuation operators 
within the different coupled channels A„k- On account of Eq. ( [4.7|) these operators take 
the form 

dr"^ r dr * 



^i-K)ij =\:r^ + -:r-ml- % - V(„K)^j. (5.4) 



The angular momenta /j, masses rrii and potentials V(^o-K)ij of these operators in the various 
channels are given in Appendix A. The fluctuation determinant decomposes accordingly. 
Therefore 



^ K=0 



det A(°) det A(°) det A^^^ 

ciet ^\jv/j^ _^ get ia^j^ ^ get i^ppx 

det' Amk det' A ek det AppK 



(5.5) 



Here it has been used that the sector consisting of the time components of the gauge field 
has the same fluctuation operator as the Fadeev-Popov sector. The minus sign in front 
of the Fadeev-Popov contribution arises from a factor (1 — 2) where the 1 correspond to 
the time components of the (real) gauge field and the —2 to the (complex) Fadeev-Popov 
ghosts. The apostrophes ' indicate that in the evaluation of the determinants the zero 
modes have to be removed and the negative squared frequency of the unstable mode has 
to be replaced by its absolute value. 



6 The Numerical Procedure 



Consider the contribution of one of the coupled channels (crK) to the logarithm of the 
fluctuation determinant n 

„ , 1 det A^^l , , 



Since our discussion will be 



with the operators A^k and A^°^ as defined in Eq. 
confined just to one channel the indices K and a are omitted in the following from the 
operators A, their Green functions, eigenfunctions and eigenvalues. 
Let us define the Green function Gij(r, r', v) by 



Aij - v^bij 



Gjk{r, r\v) 



-5{r - r')5ik 



(6.2) 
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and the analoguous Green function G^^^ for the operator A'^''^. The solution to this equation 
can be written (using discrete notation formally) as 



in terms of the orthonormalized eigenfunctions (f)ja of the fluctuation operator A. The latin 
letters refer to the field components and the greek letters label the different eigenmodes. 
Analoguous relations hold for the vacuum Green function. We define now a function F„x{^) 
by 

F.k{i^) = J drr\G,,{r,r,u)-GS\r,r,u)). (6.4) 
Using the expression (|6.3|) it can be written in terms of the eigenfrequencies Aq as 

f-M^B^-^^). (6.5) 



a 'a 



Integrating this expression over u dv from e to A yields 

/■A 1 \(0)2 , 2 ^ (0)2 , a 2 

Taking into account the fact that the operators A and A^^^ and therefore their eigenva- 
lues are linear in the ml this expression is a Pauli-Villars regulated version of (InK)^-;^, 
the second term in the sum being obtained by replacing all mf by mf + A^. Of course 
we should let A — > oo after the expressions have been renormalized and the integral has 
become ultraviolet convergent. Renormalization will be discussed in section |^. 

The lower limit e should of course be set equal to zero. We have introduced it because 
for if = 1 we have six zero modes and the limit e — ^ then obviously does not exist. We 
will deal with this problem below (see section ^) when we discuss the removal of these 
modes. 

A further problem arises from the unstable mode which leads to a pole in the region of 
integration and requires a precise definition of the integration contour. This problem will 
be tackled also in section H. 

Leaving aside these problems for the time being we have the "naive" relation 

POO 

{\nK)^K = dvvF^Kiy). (6.7) 
Jo 

In order to evaluate F^k{i^) the Green functions have to be determined. We will not use 
their expansion with respect to eigenfunctions (see Eq. (|6.3|) ), but use another standard 
method: 

Let /"^ be the solutions of the homogeneous differential equations 



'-"^ - 4) /^(r) = K„.(r)/J*(r). (6,* 
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where k„ has been defined asQ 

Kn = ^Ju^ + ml (6.9) 

whith rUn = 1 or ^ depending on the field component n. f^~^ designs the solution regular 
as r — s> oo and /"~ the solution regular as r 0. The index n labels again the field 
components, the greek letters will in the following label a set of linearly independent 
solutions of the system; there are of course as many such solutions as there are field 
components. 

For the vacuum {Vnn' = 0) these equations are solved by modified Bessel functions with 
the argument Zn = Knf, which are defined here (slightly deviating from 0] ) as 



k{z) = yll-K,^,_{z). (6.10) 



Their Wronskian is given by 

klj = Kn (il„{Zn)kl{Zn) - ki„{Zn)il{Zn)) = (6.11) 

The behaviour of the solutions for r — >■ and r — »■ cxd is analogous to the one of 
these free solutions: 

, ' 1 for z —> and ' , 1 for z — >• oo. (6.12) 

fci oc ^-('+1) J kioce-'^/z} ^ ' 

It is convenient to split off the Bessel functions from the solutions via 

n^{r) = [bl^K^{T)\h^^{z^) (6.13) 

with 

hf =H and = ki. (6.14) 

On account of the boundary conditions for the functions /^^ the functions tend 
to a constant as r — > cxo. The choice 

lim K^ir) = lim /i'"^(r) = (6.15) 

determines the normalization of the amplitudes in such a way that their Wronskian 
satisfies 

W-^r) = E (/r - fr^fnA = (6.16) 



^The letter k has been introduced previously to denote the fluctuation determinant. Since Kn as defined 
here will always appear with an index there should be no confusion. 
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For r — » the centrifugal barriers for the /"^ are different from the ones for the free 
solutions. Effectively they differ at most by two units from those, so that the functions 
5" + /i"^ behave as r°, r^^ or r^^ in this limit; most of the typical r'" behaviour that is 
dangerous for numerical calculations in high partial waves has however been taken out by 
splitting off the free solutions. 

For our new functions /i^^ equations (|6.8|) yield the inhomogeneous differential equa- 
tions 



h', izr,) d 



(6.17) 



They can be integrated numerically using the Nystrom method (Runge Kutta integra- 
tion for second order differential equations ||21[). Some numerical details are discussed in 
Appendix C. 

From these solutions the Green function defined in Eq. (|6.2| ) is now obtained as 



G^Ar. r\ v) = 0(r - r')f:^ir)Cp^f^rir') + e(r' - r)/„"-(r)C-^V„^'"'(r') (6.18) 
where the coefficients Cap are related to the Wronskian of the amplitudes /"^ as 



r'Wasir) 



>al3- 



(6.19) 



With these preliminaries we obtain for the trace of the Green functions at r' = r: 

n 

G^r^r^v) = (6.20) 



Therefore, using Eq. ( |6.13|) , the function F„kW) can be calculated for each system of 
coupled channels {crK) as 



once the functions /i"^ for that channel have been found numerically. 
The total fluctuation determinant is then obtained as 



In K 



dvvFiv) 



with 



(6.21) 



(6.22) 
(6.23) 



F{v) =Y.{2K + 1){Fek{^) + FmkW) - F^pM) 

K=0 

This expression is still formal, we have to discuss the treatment of zero and unstable 
modes and of renormalization. 
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7 Renormalization 



As has been discussed by Carson et al. [[TJ] renormalization requires, in the T — >■ cx) hmit 
discussed here, the replacement of the Higgs vacuum expectation value by its temperature 
dependent value and the subtraction of the tadpole graphs with Higgs fields as external 
legs. 



In order to do this we discuss first the relation of our expression for In k of Eqs. ( |6.21| ) 



( |6.23| ) to a Feynman graph expansion. The function F{v) is obtained as the trace of the 
Green function of the "small fields" Oj, rj and the Fadeev-Popov fields in the external 
potential generated by the classical fields. It has been obtained here by decomposing this 
Green function first into partial waves and then summing over the individual partial wave 
contributions. This Green function may also be expanded with respect to the external 
potential; formally 

°° -1 1 

G,,{x,x\v) =< x,t\y[ ^ V(x)r . . , \x',j > (7.1) 

where bold face letters indicate operators and/or matrices. Then 

F{iy) = J £xGu{c,x,p) (7.2) 

and therefore the fluctuation determinant is obtained via 

In K = j dvv F{v) 

°° \ r —1 
= E^y^'^<^^>^li^?^V(x)nf,.> (7.3) 

which is just the normal Feynman graph expansion of the 1-loop effective action for a 3 
dimensional theory. 

The tadpole contributions to F{v) are therefore obtained as that part of the first order 
contribution in the external potential that is generated by the classical Higgs field. These 
terms are easily recognized in the second order Lagrangian Cgj presented in Eq. ( |3.10| ) as 
those containing <1>^;, <1>)!^ and Phid- In the partial wave potential given explicitly in 
Appendix A these terms appear in the diagonal elements and are proportional to {H^ — 1). 

The tadpole terms can therefore be subtracted either in the single partial waves as the 
first order perturbative contribution generated be the (ifg — 1) terms or directly from F{v). 

For each partial wave the differential equation for the radial wave function may be 
transformed into an integral equation and this integral equation can be used for a pertur- 



bative expansion; this has been discussed extensively in |14[. The first order contribution 
of such an expansion is obtained as 

(ln«:)«=-/ dvvKl drr' dr'r'V„„(r')zLK<)A;f„(/s:„r;) (7.4) 
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where = min(r, r') and r'^ = max(r, r'). The tilde over the potential indicates that only 
the terms proportional to {Hq — 1) are to be included. 

One of the radial integrations can be performed analytically so that 

\ /"oo poo _ 

= - / duu rfrrV„„(r) [ii^-ik^ - k^+ikj ■ (7-5) 
2 Jo Jo 

Either by summing up the partial wave contributions or by calculation of the Feynman 
type graphs according to Eq. {TA_) the total tadpole contribution to -F(z^) is obtained as 

FU^) = -- drr\H',{r)~l)[^ + ^ + ^] (7.6) 
2 Jo Kw 2kh 2kw 

where Kn has been defined in Eq. ( |6.9D . The renormalized value of In k is therefore obtained 
from 

= F(z/) - F,,,(z/) (7.7) 

where the subtraction can be done either in the partial waves using Eq. (|7.5| ) or in the full 
amplitude using Eq. (fTOl). The z/ integration of Freniy) is now ultraviolet convergent, i.e. 
the upper limit of integration A introduced in Eq. ( p.6|) which serves also as a Pauli-Villars 
regulator can be sent to oo. 



8 Zero and Unstable Modes 

Besides the question of renormalization we have also postponed the treatment of the zero 
and unstable modes. 

The zero modes should be removed from the fluctuation determinant. Equation ( |6.(j| ) 
shows that their contribution to {\yik)ei and (Infi;)^! is in both cases — (l/2)ln(AQ + e^) 
before e goes to 0. Of course Aq = here and the limit does not exist. But we have to 
remove just these contributions. Therefore the fluctuation determinant without the zero 
modes is obtained as _ 

poo 

lnK = lim(/ rfz/z/F^e„(z/) +61n(e)). (8.1) 

The function Fren{i^) will of course behave as due to the zero modes so that the 
limit exists. It has to do so also in the numerical evaluation; this represents a good cross 
check. Of course e is a quantity of dimension energy and indeed removing the zero modes 
makes k a quantity of dimension {energy^. This has to be taken into account when 
comparing results if different length and energy units are used. We have used units of 
for the radial variable; our unit for the eigenvalues of the second order differential 
operators A is therefore M^r ,their eigenfrequencies and therefore also the variables v and 
e are therefore in units Mw- 

The unstable mode is not to be removed but, according to the general theory ||, 0] to 
be replaced by its absolute value. If its eigenvalue is denoted by = —ui"^ < then it 
leads to a pole in F{v) at z/^ = tu^ > and therefore in the region of integration of Eq. 
( |6.22| ). It would contribute a term — (1/2) ln(— cj^ ) and the minus sign in the logarithm 
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should be removed. This can be done simply by evaluating the integral as a principal value 
integral. 

Since doing singular principal value integrals numerically is a delicate operation we 
have subtracted the pole from the integrand and done its integration analytically. We have 
used the identity 



(Iuk) 



P / duvF{v) 
Jo 



dv V 



— up_ 



+ 



Z/2 + 



— ln(u;_/cr). 



(8.2) 



The third term in the parenthesis was added in order not to spoil the ultraviolet convergence 
of the integral. Since (^.2|) is an identity the value of a is in principle arbitrary. We have 
chosen cr = 1 in units of Myi/ for convenience. 

For easier notation have treated the removal of the zero modes and of the unstable 
mode singularity separately. It is of course understood that both prescriptions Eqs. ( |8.1| ) 
and (|8.2|) are applied simultaneously. 



9 Results 

We have now discussed the principles of our numerical procedure for calculating the fluc- 
tuation determinant. Before presenting the results we will discuss some specific details of 
our numerical evaluation. 

The first step was the Runge-Kutta integration of the partial wave differential equations 
( |6.1?1 ) in each channel and the numerical integration of the exact trace as presented in Eqs. 
( |6.21| )-( |5.2^ ). The numerical details are discussed in Appendix C. 



In order to obtain the function F{v) we had to sum over all partial waves K within the 
various sectors a and then to perform the summation over a. In order to have a check on 
the K summation we have considered the asymptotic behaviour at large K of the terms in 
this sum. For this purpose it is sufficient to consider the perturbative contribution of first 
order in the potentials V^j, higher order contributions will decrease faster. The first order 
contributions to the sum in Eq ( |6.23|) have the form 



{2K + l)Fi]} = {2K + 1) / dr-VUr) K-iK. - ■ (9.1) 



Using the uniform asymptotic expansion of the Bessel functions at large K one finds 
that this expression behaves as asymptotically. This determines then the convergence 
behaviour of the K summation. It has been checked numerically to a good accuracy. Since 
we know the leading behaviour, we can extrapolate the terms to arbitrary valus of K. We 
have done so by fitting the last five calculated values with a power behaviour 
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To the sum extended to some value K^ax 

F{K^ax, i^) = E E(2^ + (9-3) 

K=l cr 

we have then added the sum from {Kmax + 1) to oo using the fit (|9.2|). Fig. 3 shows 
the partial sums as well as the sums completed using the extrapolated values. One sees 
that the complete sum becomes independent of Kmax already at moderate values of this 
variable. This shows that the extrapolation procedure is reliable. 

The function -F(z/) is displayed in Fig. 4 for u = 1. The pole contribution of the unstable 
mode is already subtracted here according to Eq. (p. 21). The dashed line shows the full 
function F{i>), the dash-dotted line the tadpole contribution. Note that this contribution 
is determined analytically in absolute normalization. So the fact that both curves approach 
each other as u ^ oo checks also the absolute normalization of the unsubtracted F{u). 
The full line shows Freni^^) together with the asymptotic estimates at small and large u 
(dotted lines). The behaviour at small u is normalized absolutely; it is determined by the 
zero modes to be 

The u integration was performed numerically up to Umax — 2.5, then an asymptotic 
part was added to the integral by extrapolating Fren as Ci/ + Di/ k^. The results for In k 
are given in Table ^for various values of ^ = Mh/Mw in the scale Mw- They are plotted 
in Fig. 5 together with previous results and estimates in the scale gv, i.e. after subtracting 
6 In 2 from the values of Table |l|. Our calculation stops at Mh/Mw = 2 for a technical 
reason: Above this value the leading asymptotic behaviour exp(— K^^r) of the fluctuation 
hi of the Higgs field becomes dominated - through a cross term in the potential - by a 
gauge field contribution which behaves as exp(— (M^y + kvk)^)- So for M^ > 2M]y the 
boundary conditions for this function have to be modified. If the Higgs mass should turn 
out to be larger than this value we would have to deal with this complication. 

We think that the main uncertainty of our results comes from the extrapolation of 
Fren{i^) to u = OO. Wc estimate the error of this asymptotic contribution to be around 
10% yielding an typical error of 0.3 for In k, and we think that this is a conservative estimate. 



Our results differ considerably from the ones of CLMW [0]. There is one point which 
could lead to a difference on physical grounds: we have used a different gauge in the classical 
sphaleron ansatz (see below Eq. ( |2.17D ). While the sphaleron transition rate must be gauge 
invariant, a difference in the classical configuration could modify the zero mode prefactors 
and then require also a compensating change in the fluctuation determinant However 
- as mentioned above - the prefactors agree and therefore should also the fluctuation 
determinants. 

The methods used by CLMW and by us differ also considerably. Since we work with 
the Euclidean Green function, we have no difficulties with the fact that the spectrum of 
the various fluctuation operators is continuous (except for zero and unstable modes and a 
few further bound states). The Schwinger proper time method as used by CLMW requires 
a discrete spectrum which has to be created artificially by introducing a space boundary 



^We thank L. McLerran for pointing this out to us. 
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(a kind of "bag") of radius R. The calculation requires then a limit R ^ oo. This limiting 
procedure is absent in our method. This makes the algorithm much faster. Furthermore 
the method contains more internal consistency checks. Of course this does not necessarily 
imply that our results are the correct ones but at present we have not found any reason to 
doubt them. We hope that the discrepancy can be settled in the near future, possibly by 
another calculation. 

Comparing to the DPY approximation we find that our results come, at Mh > My/, 
much closer to this approximation than the ones of CLMW; they show a different trend for 
small Higgs masses, however. The validity of the DPY approximation has been discussed 
by Carson for a one-dimensional sphaleron where the exact fluctuation determinant 
is known analytically. He finds only fair agreement between the exact results and the 



approximation. Furthermore (see |]22|) it is not obvious in which way the convergence of 
gradient expansions (and the DPY approximation falls under this category) in 1 + 1 and 
3 + 1 dimensions can be compared. 

In conclusion we have presented here a new set of fluctuation equations for the the 
electroweak sphaleron at Qw = 0. We have used it to evaluate the fluctuation determinant. 
We obtain the result that the quantum corrections lead to an enhancement of the sphaleron 
transition with respect to the estimate obtained from the classical saddle point solution. 
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A The Partial Wave Fluctuation Equations 



The equations of motion have been obtained as follows: the Ansatz for the fluctuation 
amplitudes was inserted into the second order fluctuation Lagrangean which itself had been 
obtained using REDUCE from the general gauge-fixed Lagrangean. A second REDUCE 
step used the various algebraic properties of the operators I and J given in in order 



to obtain a Lagrangean quadratic in the amplitudes ta, Sa, Pa and h. The Lagrangean 
before gauge fixing had already be obtained in ||l3l . Checking its gauge invariance and the 
presence of all zero modes presents a good test on its correctness. 

The equations of motion obtained in this way are - for K > 1 - given by 

rH'l + 2rt[ - r% 

= hir^H^ + 2/1 + /C^ + 2) - 2t2fAlC' + 2h{2rf'^ - 2/^ - /C^) 

- 2UlC\rfj, - /a - 1) - 2t,ie{rfj, - /a + 1) + 2p,r^H', 
r% + 2rt'^-rH2 

= - 2tjA + Ur^Hl + fl + K? + l)- 2h{rf\ - /a - 1) 

+ 2t6(r/; - /a - /C^ + 1) - 2t,{rt\ - /a + 1) + 2p2r^H'^ 

rhl + 2rt'^-r% 

= ~2t,- 2t2{rf^ - /a + 1) + hir'H'o + 2/1 + /C^) 

- t6[(3/A - 1)(/a - 1) + 2/a/C2] + t9(3/A - 1)(/a - 1) 

- P2rHo{fA - 1) 



rHa + 2rt'-rH 



6 



= {-2ti(r/l-/A + l)-2t2(r/l-/A + /C2 + l) 

- h{fl + 2/a/C2 - 1) + h[2iefA{fA - 1) + (3/a - 1)(/a - 1)] 
+ h[ie{r'Hl + /I + - 1) - (3/^ - - 1)] 

- pirH^{fA - 1) - P2rHo{fA - 1)}//C' 



rhl + 2rtg-r% 



= {- 2ti(r/l - /a + 1) - 2t2(r/l - /a + 1) - h{fl - 1) 

+ t6[2;CVA(/A-l)-(3/A-l)(/A-l)] 

+ t,[]C'ir'H', + /I + _ 1) + (3/^ _ _ 1)] 

- p,rHo{fA - 1) - P2rHo{fA - 1)}/IC' 



r 



p^ + 2rp^ — r pi 



= {ihr^H', + 4hrHo{fA - 1) - 2terHo}C\fA - 1) 

- 2t,rHo}C^ifA-l)-2p2lC\fA + l) 

+ pi[r^e{Hl - 1) + 2r^Hl + (/a + 1)^ + 2ie]}/2 

r'^P2 + 2rp2 — r^p2 

= {At^r^H'^ - 2tsrHo{fA - 1) + 2UrH,{fA - 1) 

- 2^rH^{fA-l)-2pi{fA + l) 



19 



+ p2[rY{H', - 1) + 2r'H', + /I + 2/0^ - l]}/2 (A.l) 

for the magnetic sector, by 

+ 2t7(r A - - /C^ + 1) + 2t^{rfj, - /a + 1) + 2p3r'i/^ 
rH'i + 2rt'5-r2t5 

= -2t4(r/;i-/^ + l)+t5(rX + 2/l + /C2) 

- t7(3/A - - 1) - tsP/A - - 1) + 2/a/C2] 

- p3rHo{fA - 1) 
r^tl,' + 2rt'^-r% 

= {2U{rf^-fA-ie + l)+h{fl-l) 

+ h[ie{T^Hl + /2 + /C2 _ 1) + (3/^ _ _ 1)] 

- t8[2/CVA(/A-l)-(3/A-l)(/A-l)] 
+ pzrHo{fA - 1) - /iiri/o(/A - 1)}//C' 

rHl + 2rt'^-r% 

= {-2t4(r/;i-/A + l)-t5(/l + 2/A/C2-l) 

- i,[2X:VA(/A-l) + (3^-l)(/A-l)] 

+ t8l/C2(r^i^o' + /I + /C^ - 1) - (3/a - - 1)] 

- P3rHoifA-l) + hrHo{fA-l)}/K:'' 
r^Ps + 2rp3 - r% 

= {4i4r'//o - '^hrHoifA - 1) + 2t,rH^{fA - 1) + 2tsrHo{fA - 1) 

+ P3[r'e'(^o - 1) + 2r^/^o' + /I + 2/C^ - 1] 
+ 2h{fA-l)}/2 

r^h'l + 2rh[-r% 

= {- At.rHoifA - 1) - 2tjrHo)C\fA - 1) + 2tsrH^lC\fA - 1) 

+ 2p3/C'(/A - 1) + /iiPf (3//o - 1) + (/a - If + 2/C']}/2 (A.2) 

for the electric sector and by 

2 II I o ' 2-- 

r + zrS]^ — r Si 

= si (r^ifo' + 2/1 + X:^) - 2s2/a/C2 

r S2 + 2rs2 — S2 

= -2si/^ + S2(r2i/2 + /2+K:2-l) 

r^Sg + 2rs3 - r^ss 

= s3{T^Hl + fl + K?-l) (A.3) 

for the Fadeev-Popov sector. 
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As one sees these equations do not have the form of Eq. ( |4.7| ) with a symmetric potential 
Vij. Such a form is however needed on order to apply the Green function formalism in a 
convenient way. 

The transformations that bring the equations of motion to a symmetric form have been 
found by an - in fact very simple - educated guessing. They are 



K{K+l)r {K + 2){K+r 



V2K + T 



2K + 3 



K 



V2K-1 



To 



lK{K 



K 



V2K + 3 



2K 



+ 



{2K + l){2K + 3) 
K+1 



{K+l)Ts-J{K+l){K + 2)T, 



{2K~1){2K+1) 
K{K+1) 



{2K-1){2K+1) 

K{K + 1) 
(2K + l)(2K + 3) 
K+1 



K{K - 1)T6 - KT2 
T3 + 



/ K 
K-1 

lK+1 



K + 2 



{2K + 1){2K-1) 
K 

{2K + l){2K + ?>) 



T2 + 
T3 + 



K 



K-1 

lK+1. 
K + 2 



-n 



Ti 



^9 
Pi 

P2 



K{K + 1) 
V2K+1 



K{K+1) 
V 2K + 1 

for the magnetic amplitudes and 

1 

h 

tr 
U 



VKP2 - VkTTPi 

VkPi + VKTIP2 



V2K+1 
1 



1 



VkTT 



K 



V2K~+T^V1<^ Vl<T2 



K{K + l)^/2K +1 



P3 

hi 



V2K + T^{K + l)VK + 2 

^ Ps 
K{K+1) 

Hi 



Kn + y/KTTn 



Ky/lT 



(A.4) 



(A.5) 
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for the electric amplitudes. For the Fadeev-Popov amplitudes the transformations are 

(A.6) 



Sl 

S3 



^fKS2 - yFTT^i 

5*1 ;^»S'2 



S3- 



The n-tuples formed by the amplitudes of the various sectors 



FP\ 



(Tl,T2, T3, Tg, Tg, Pi, P2] 

iT,,T,,T-r,Ts,P3,H,) 
{Sl, S2, S3) 



(A.7) 

- the 



fulfill - with the mass parameters rrii and the angular momenta /j from table | 
required fluctuation equations ([4.7| ). 

The symmetric potential of the fluctuation equations in the form (17) has for the 
magnetic sector the elements 



yM 



+ 



fA 



V, 



M 
22 



{Hi - 1) + 



M 
33 



4/A(^ + 2) 
r(2ir + 3) 

4ir2 + 3(/^ + 5)ir + 7/A + 15 

^fA{K - 1) 
r{2K-l) 

8K^ - Q{fA - 1)K' + {5fA - 7)K^ - AK ~ 3fA + 1 

/a-1 



1 



X 



yM 



r{2K + 3) r^{2K + 3){2K +1){K +1) 
8K^ + 2(3/a + 13)^3 + 23(/a + l)K^ + 4(7/^ + 1)K + 8/a 
4/^(^-1) 



(Hi - 1 



M 
55 

M 
66 



r'^{2K-l 
{Hi - 1) + 
1 



r{2K -I) 
fA-l 



2\/u2 



[{fA + l)K{K + l) + 3fA-l] 



^,2 + e){Hl-l) + ^^{2K + fA + 3) 



77 

M 



i(2+a(i^o^-i)-^ 



1 



(2ir - /a - 1) 
iK+l)iK + 2)i2K-l) ifA-lf 



2K + 3 



r^{2K+l) 



22 



^ K+l [r{2K + 3) 
r^{2K + 3){2K +1) 



[4{fA + l)K + fA + 3] 







K + 2 



X 



{2K + 1){2K + 3)(K + 1) 
2f^K U - 1 



[2{fA+l)K-fA + l] 



^^^;-[2rH',-(U-l)Ho] 



2K + 3r 



2{fA-lf 



^{2K + 3){2K -l)r\2K + 1) 

iK-i r 2f'^ u - 1 



K [r{2K-l) r2(4X2_i-) 
1 {2f^{K+l) 



[A{fA + l)K + 3fA + 1] 



V4.K'^ - 1 1 



+ 



2(/^ - 1)^2 + (/^- 3)^-3/^ + 1 



/ir(A--l)(2A' + 3) (,f_i - 1)2 



2is:-l r2(2is: + l) 



^(2X + l)(2X + 3) 
2/^X , - 1 



X 



+ 



2{fA - l)K^ + {3fA - 1)K -2fA + 2 



K{AK'^ - 1) 



23 



v,^ = 



yM _ 
^47 — 



yM 



I K-1 1 
2K-lr 
K 



[2rH', - if A - l)Ho] 



M 
57 



^{K+l)(2K+l)r 



^K{2K + 1) ^ 



M 
67 



0. 



For the electric sector we obtain 



- (^0-1) 



yE 



33 



44 



1^. 



55 



'^66 



yE 
^12 



Vio — 



13 



14 



^1 



15 



r(2/s:+ 1) 



r2(2i^ + 1)(X + 1) 
2 .^ , 4/A 



2(/^ + l)K^ + (3/^ + 7)ir + 4/^ + 4 



r(2/s: + 1) 



(^0^ - 1) + 

r'^K{2K^l) 
{Hi - 1) 

+ 

r2(2is:+ + 1) 

-2 .^ 



(^0 - 1) 



2{fA + + (/a - m + 3/a - 1 

4ir3-(3/A-l)(K^-X-l) 



r'^K{2K+l) 



l(e^ + 2)(//o^-l) + ^(/l-l) 



1 



2r/M-(/^-l) (3/^+1) 



^K{K + 1) r2(2/r+l) 
^K{K + 2) (2f^ , /^-l 



2X + 1 



r ^r2(X+l) 
-{2r/;i(X+l) 



[2(/^ - 1)K -fA-l] 



K r2(2K + l) 
-(/^-l)[2(/^ + l)X + (3/^ + 1)]} 



1 



^{K + l){2K + l)r 



-[2rH',{K+l)-{fA-l)Ho] 



24 



yE 



Vr 



/ K 1 
IK + 2 



if A - l)Ho 



23 



K +lr^{2K +1) 



{2rf'^K - if A - 1) [2ifA + l)K -fA+ 1]} 



24 



is:(2is: + l)r2 



{2r/;iK + (/^ - 1) [2{fA - l)K + 3/^ - 1]} 



25 



Vr 



26 



•^34 



2/^ + Ir 



^{K-l){K + 2) 
{2K+l)r^ 



if A - If 



35 



yE 



yE 



^ (X + l)(2X + l)r 



2K + lr 



-1 1 

-(^-l)i/o 



^ X(2K + l)r 



46 



//r-i 1 

2K + lr 



if A - l)Ho 



yE _ 
^56 — 



(A.9) 



and finally for the Fadeev-Popov sector 



FP 
11 



yFP _ 



yFP 
'^33 



(^0 - 1) 
(^o' - 1) 



fA-l 

{2K + iy 

fA-i 

{2K + l)r2 



47^2 + (3/^ + 7);^ + 2/^ + 2 
4/^2 - (3/a - l)K -fA-l 



(^o-l) + -^(/l-l) 



12 



yFP 

FP 



{2K + l)r2 



13 



FP 
23 







0. 



(A.IO) 

We have again used REDUCE to obtain these potentials from the equations of motion 
for the transformed amplitudes. 



25 



In the case K = the I- and J-operators act on the constant spherical harmonic loo- 
Therefore only the operators I^, I^, and contribute to the fluctuation Lagrangean. 
We can account for this by setting - in the Lagrangean - the amplitudes t2, t^, t^, tj, tg, 
tg, p2, ps, S2 and S3 equal to zero. For the remaining amplitudes ti, t^, t^, pi, hi and si 
one obtains then the Euler-Lagrange-equations 

rh" + 2rt[ - rHi 

= ti{r^Hl + 2/1 + 2) + 4t3(r/; - /^) + 2pir^H', 

rHl + 2rt'^-r% 

= 2ti{rf^ - Ia) + h{r^Hl + 3/1 - 1) + pirH,{fA - 1) 

rH'i + 2rt'^-r% 

= hir^H^ + 3fl~l)-hirHoifA~l) 
r'^Pi + 2rp[ — r^pi 

= {4tir^H', + At,rHo{fA-l) 

+ PiyeiHl - 1) + 2r'Hl + {fA + l)^]}/2 

r'^h'l + 2rh[ - r'^hi 

= {-At.rHoifA - 1) + hilr'ei^H', - 1) + (/^ - lf]}/2 

2 II I n I 2- 

r Si + 2rsi — r Si 

= si{T'Hl + 2fA). (A.ll) 



Changing the basis as 

ti = V2T3 - 2Ti t, = V2T5 si = Si 

ts = V2T^ + Ti hi = 2Hi (A. 12) 

Pi = VQPi 

and introducing the n-tuples 

(r/f) = (ri,T3,Pi), {vf) = iT,,Hi) and r/^^ = (A.13) 



for the fields the differential equations take the required form of Eq. The masses rrii 
and angular momenta /j are shown in table ^ 
The potentials in the three sectors are 

8r/A-(/A-l)(7/A + 15)] 



1 

2? 



yM _ 
Vii — 




yM _ 
^22 — 




yM _ 




yM _ 
^12 — 





■^-l)(/^ + 3)] 
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V^^ = -^^[{fA-l)Ho + rH^] 

Vlf = - — ifA-l)Ho 
r 

V^'' = {Hi - I) + -Ml- I). (A.14) 



Another exceptional case arises for = 1. Here the J- and I-operators act on the 
spherical harmonics Yim x. The action of becomes equal up to a sign to the one of 
and similarly the one of to the one of (see Eq. ( |4.3| )). Therefore a linear combination 
of the amplitudes ts and and of ts and tg respectively can be chosen to vanish. We have 
taken the choice 

= and tg = 0. (A.15) 
For the remaining 14 amplitudes we find 

rh'l + 2rt[ - rHi 

= h{r^Hl + 2/1 + 4) - AhfA + 4t3(r/^ - /a - 1) 

- At,{rf^- fA + l)-2p^r^H'^ 
rHl + 2rt'^ - rH2 

= -2tjA + h{r^Hl + /I + 3) - 2h{rfA -fA-l) 

- 2t,{rf^- fA + l) + 2p2r^H'^ 
r% + 2rt'^ - r% 

= {2t,{rf'^- fA-l)- 2h{rf'^- fA-l) 

+ h{2r^Hl + 5/1 + AfA + 3) - t,{fl - 1) +pirHo{fA - 1) 

- P2rHo{fA-l)}/2 
rHg + 2rt'g-r% 

= {-2ti(r/l - /a + 1) - 2t2(r/l - /a + 1) - Wl - 1) 

+ t,{2r'Hl + 5fl-AfA + 3)-p,rHo{fA-l) 

- P2rHoifA-l)}/2 
r^p'l + 2rp[ — r'^pi 

= {At.r^H', + AhrHoifA - 1) - At,rH^{fA - 1) 
+ pi[r^e{Hl - 1) + 2r^Hl + /I + 2/^ + 5] - Ap^ifA + l)}/2 
r'^P2 + 2rp2 — 'r'^V2 
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= {2i2r'//o - '^hrHoifA - 1) - 2UrHo{fA - 1) - 2pi(/^ + 1) 

+ P2[r'e{H', - 1) + 2r'H', + fl + 3]}/2 (A.16) 



for the magnetic sector, 

rH'l + 2rf' - rH: 



4 



= Uir^H^ + /I + 3) - 2t,{rf'A - /a + 1) + 2i7(r/; - /a - 1) 

rHl + 2rt5 - r% 

= {-2U{rf'^ - /a + 1) + t5(2r2if2 + 5/2 _^ ^ 3) 

+ trifl - 1) -P3rHo{fA - 1) - /iiri/o(/A - l)}/2 

= {2U{rf'^ -fA-l)+ Wl - 1) + tr{2r^H^ + 5/1 - 4/^ + 3) 

+ psrHoifA - 1) - /iiri/o(/A - l)}/2 

r ^3 + 2rp3 — r p3 

= {UyH'^ - 2t5rHo{fA - 1) + 2trrHo{fA - 1) 

+ P3[rT(i/o' - 1) + 2r^iil + /I + 3] + 2/ii(/a - l)}/2 

r^/ii + 2r/i; - r^/^i 

= {-4t5ri/o(/A - 1) - 4t7ri/o(/A - 1) + ^Vz{iA - 1) 

+ /iipe2(3i/o'-l) + /l-2/A + 5]}/2 (A.17) 

for the electric sector and 

r^s'[ + 2rs\ - r^Sy 

= si(rX2 + 2/l + 2)-4s2/A 

r^S2 + 2rs2-r^S2 

= -2si/A + S2(r2i/2 + /2^1) 

2 // I o ' 2 •■ 

r S3 + 2rs3 — r S3 

= S3(rX + /l + l) (A.18) 

for the Fadeev-Popov sector. The amplitudes are transformed as 

= (2T1 + V5T2 + Vera) ^4 = ^ (V2T5 - 2r4) si = 52-V25i 

^2 = (Vers + 2a/5T2 - 2Ti) = S2 = -S2 + -^-Si 

^3 = (v^Ta + STi) ^ _L + ^Ts) S3 = 53 (A-19) 

U = VSTg P3 = V2P3 

pi = 2 (P2 - V2A) hi = 2//1 

P2 = 2P2 + V2P1 
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in order to get a symmetric potential and asymptotic decoupling. With the n-tuples 



2) 



ivD = iSi,S2,Ss) (A.20) 

the fluctuation equations take the form of Eq. ( [4.7| ). The masses rrii and angular momenta 
li are given in table 

The elements of the symmetric potential are given by 

rM i„i ,^ , , (//i-l)(65/^ + 61) 



T/M _ (tt2 in, ^Ua — 1) 
^22 - (^0 - 1) + ^2 

.rM ,rr2 1 2/;, (/^ - 1 ) ( 10/^ + 34) 



5 

(/a-1)(5/a + 1) 
2r2 



2(/a-1)^ 



•^12 

1 



3V5r 



V,^ = ^^[(/^-l)(5/^ + 7)-6r/^ 



5-\/6r 



AJ 



= --^[(/.-l)(3/.-l) + 4r/;] 



^if = -^mA-l)Ho + ArHl,] 
vlOr 



= 



yM 



2 (/a - 1) 



2 



23 \/ ^5 ^2 



= ^[(/A-i)-r/;,] 



= 

^26^ = 2i7o 

rM 1 



^34 = --^[(/A-l)(/A + 3)-2r/;,] 



Vi^ = -Jll[ifA-l)Ho-2rH',] 
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— n 

1 



V,^ = -^(^-l)//o 

= (A.21) 
for the magnetic sector, by 

rrr2 ^^ "^Ja , (/a - 1) (9/a + 13) 

(rr2 1^^4/;, (/^-l)(6^-2) 
^^22 = (-^0 - IJ + ^ + 



^3! = (H^-l) 



3r Sr^ 

(./.i-l)(5./,i + 9) 
2r2 



^/E ■'■/'O I t2\/'o-2 I (/-A + 1) 

^44 = 2^ ^ ° ~ ^ 2^^2 

^55 = 2^ ^ ~ ^ ^ 2^^2 

V,^ = _^^[2r/;-(/^-l)(3^ + l)] 



3V^r2 

nf = ^^[4r/U(/A-l)(/A-3)] 

= --i-(/A-i)i/o 

^2! = -^[2r/;i-(/^-l)(/^ + 3)] 



^2! = [2ri/^ + (/a - l)i^o] 

^3! = --(/a-1)//o 



= (A.22) 



for the electric sector and by 



V^P (H^ , (/a^:1)(5/a + 13) 
30 



V- = K - 1) + :*(-^^ - 



rFP 

13 

tFP 



3r^ 



Vi/ = (A.23) 
for the Fadeev-Popov amplitudes. 
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B Zero Mode Prefactors 



The prefactors Mtr and Mrot are determined by the normahzation of the translation and 
rotation zero modes. We have included factors l/27r which otherwise would appear with 
the prefactor T^^. We have then (cf. 0) 

Mtr = Nl 

KotVrot = StT^NI, (B.l) 

where the normalization factors are given by 

^ J dW:'^°Y:'^°'- (B.2) 

The ipn are the zero mode wave functions, the different field components are assumed to 
have canonical normalization (i.e. appearing as ^d^ipj^d^ipn for each real component in the 
Lagrangean density). The rotation and translation modes have been determined explicitly 
T3[| for the sphaleron solution in the form ( p.l7| ). The zero mode amplitudes were found 



m 



to be 



fA + l f'A 



A 



7 



J.2 2f 
-f'A 



for the translation mode and 



2r 

/ii = H'^ (B.3) 



" ^ (B.4) 



r 

for the rotation mode. 

For proper normalization (which was irrelevant in |13[) all these contributions have to 
be multiplied by y^47r/3 which comes from the different normalization of the Yim and the 
xm used in that calculation. Also the translation mode amplitudes have to be multiplied 
by Mw if they are generated by the ordinary gradient, i.e. the derivatives w. r. t. x^. We 
will include these additional factors at the end (see Eqs. ( [B.ll| ) and ( |B.12| )). 

In this form the modes are not normalizable and do not satisfy the background gauge 
condition ( |3.7| ) . The general form of the infinitesimal gauge transformations (which also ap- 
plies for finite ones since we have expanded the fields linearly around the classical solution) 



has also been given in . For the K = 1 channel it reads: 

= g'l 

6t2 = 92 
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Sp2 

Shi 



gi + 92 if A - 1) 



2r 



(^1 + 3^2) 



93 



2 



93 



+ 



r 

fA 



9s 
1 



2r 

Hq92/'2 
Ho93/'^ 
K 



2r 

-(91+92) 



(B.5) 



The background gauge conditions read in terms of the partial wave amphtudes 

rt[ + 2ti = 2{fA + l)h-2{fA-l)t9 + 2rHoPi 

rt'^ + 2t2 = -{fA + l)t3- {fA-l)t9 + 2rHoP2 

rt'^ + 2U = -{fA-l)h + {fA + l)t7 + 2rHoP3. (B.6) 

Inserting the amphtudes and gauge functions leads to three differential equations for 
the gauge functions 



92 + -92 
r 



93 + -93 



Hi 



2(/l + l) 



4/a 2{fA 
9i - —^2 - - 



1)^ 



/l + l 



92 



93 - 



r 

Va 
~^9i 

if A - If 



(/.I 



1? 



(B.7) 



which have to be solved with boundary conditions that ensure the normalizability of the 
zero modes. The solution for 513 can be found explicitly: 

l-fA 



93 = 



(B.8) 



so that the translation mode becomes 

h -- 

h -- 

P3 = 

hi -- 



f'A 



r 

/I 



2r 
2r 



2r2 
fA - l)Ho 



2r 



(B.9) 
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The gauge transformed rotation mode becomes 



*i = g'l 
h = 92 



+ -1^ ^1 - 92} 

r 2r 

(^1 + 92) + 



2r 

Pi = 2 ^0^1 

P2 = \n^92. (B.IO) 

We have not been able, however, to solve the equations for the gauge functions g\ and g2 
analytically. The normalization integrals A^**"'^"* are obtained by inserting these amplitudes 
into the general expression amplitudes as 



Here the first factor in front of the integral is a factor "borrowed" from the factors \/l'nT 
which arise for each extracted zero mode and which we have included here as it was done in 
The factor 47r/3 has been explained above. Finally the factors in front of the different 
amplitudes come additionally from the normalization of the tensors used in the expansion. 
Explicitly we obtain 



8M, 



w 



tr 



dr 



2 Tji'i , 1 TT2t f- 1 \2 



rot 



35-2 Jo 



2 J 1 „/2 , „/2 



9i- 92 . Ia - If X 

+ -1^ [91 - 92) 

r 2r 



+2 



fA-l 

2r 



{9i + 92) 



The expression obtained for the translation mode can be shown to be 



(B.12) 



(B.13) 



as expected from a general virial theorem. The rotation mode normalization should be 
related, by a similar virial theorem, to the moment of inertia of the sphaleron. We have 
evaluated the normalization integrals numerically; we find after taking into account factors 
2±3/2 different scales M-^ and gv used in the two publications - the same results 

as within the numerical accuracy, though we have used a different gauge for the classical 
solution. Note that the scale factors cancel in the product Mrot^trans- 
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C Determination of the Amplitudes h 



n 



In this Appendix we will discuss briefly the numerical evaluation of the functions /i"^, i.e. 
the solutions of the differential equations (|6.17|) . 



The numerical integration of the differential equations ( |6.17|) is started, for /i^^, at 
some sufficiently high r = Vmax with the initial condition h'^'^irmax) = 0. 

Starting the functions at r = is by far more difficult: the behaviour of these 
functions as r — > has to be determined analytically. This means that all the functions 
that enter the differential equations (|6.17|) , i.e. the Bessel functions and their derivatives, 
the classical profiles /a and Hq and the solutions h^" have to be expanded, for r — 
into Taylor series. As a first step which proves to be nontrivial one has to find the leading 
behaviour of the solutions, since the centrifugal barrier at r = differs from the vacuum 
sector one. If we write the leading behaviour as we find 



A 



2 for the amplitude 5*3 

1 for P2 

forT5,T6,T8,52 (C.l) 

-1 f0TPuP3,H, 

-2 for Ti,T2,T3,T4,T7,T9,5'i. 



With these parameters one obtains in each n x n sector a set of recursion relations for the 
next-to-leading coefficients and a set of starting conditions for n independent solutions 
labelled by a = 1, n. 

This expansion which has been determined up to the second nonleading order in is 
used up to some suitable r = r^m at which the Nystrom integration is then started. The 
solutions found in this way do not yet satisfy the boundary condition h'^~ — >■ as r 00. 
However a set of such solutions can now be found by a simple linear transformation. 

A good check on the accuracy of the numerical integration consists in checking the 
constancy of the product r'^W°'^{r) related to the Wronskian (see Eq. ( |6.16|) ). For 

J 5 ■ 10~^ for the amplitude 5*3 ,„ . 

1 0.1 for all other sectors 

this expression was found to be constant to 5 significant digits. For smaller r the numerical 
integration becomes delicate for all sectors except the amplitude S3, since some of the 
amplitudes become singular as r — »• 0. In this region we used the known leading behaviour 
of the exact Green function 

2^ / N I for the amplitude 6*3 ,„ „^ 

I r for all other sectors 

with coefficients determined from the numerical results in the reliable region (|C.2|) . 
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Figure Captions 



Fig. 1 The classical splialeron energy Ed'. 

The figure shows the classical sphaleron energy Ed as a function of ^ = Mh/Mw in units 
Fig. 2 The zero-mode normalization factors: 

The solid line shows the normalization factor A/tr of the translation mode and the dashed 
line the normalization factor A^-ot of the rotation mode as a function of ^ = Mh/Mw- The 
units are {Mw/g'^Y^'^ and {Mwq'^Y^^'^ respectively (see Eqs.( p.l2|) ). 



Fig. 3 The convergence of the K summation: 

We show the partial sums F{Kmax, ^) as defined in the text a function of Kmax for different 
values of v. The dashed lines are the values obtained by including the sum from {Kmax + 1) 
to oo using the fit of Eq. (|9.2| ). These values are seen to become independent of K^ax 
already around K^ax ~ 10. 

Fig. 4 The function F{v) for i = Mh/Mw = 1: 

The solid line shows uFj-eni^), the dashed line the unrenormalized vF{v). The pole contri- 
bution of the unstable mode has been removed (see Eq. (^.2D). The dotted lines show the 
expected power behaviours at small and large v and the dash-dotted line the analytically 
known tadpole contribution to vF{v). 

Fig. 5 The fiuctuation determinant: 

The circles are our results, the crosses those of CLMW. The full line is the estimate of 
Carson and McLerran based on the DPY approximation and the dashed line a perturbative 
estimate. 
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Tables 





0.4 


0.5 


0.6 


0.8 


1.0 


1.5 


2.0 




1.32 


1.36 


1.39 


1.45 


1.51 


1.62 


1.71 


In K 


6.18 


5.89 


5.68 


5.50 


5.48 


5.62 


5.71 



Table 1: The results for In k for various values of ^ = Mh/Mw together with the frequen- 
cies a;_ of the unstable mode 



M 


ITii 




1 


1 


K + 2 


2 


1 


K 


3 


1 


K 


4 


1 


K-2 


5 


1 


K 


6 


1 


K + l 


7 


1 


K-1 



E 


nii 




1 


1 


K + l 


2 


1 


K-1 


3 


1 


K + l 


4 


1 


K-1 


5 


1 


K 


6 




K 



FP 


m. 


k 


1 


1 


K + l 


2 


1 


K - 1 


3 


1 


K 



Table 2: Masses and Angular Momenta of the Amplitudes for K > 1 



M 




k 


1 


1 


2 


2 


1 





3 


1 


1 



E 


nii 


k 


1 


1 


1 


2 








FP 


rrii 


h 




1 


1 



Table 3: Masses and Angular Momenta of the Amplitudes for = 
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M 


nii 




1 


1 


1 


2 


1 


1 


3 


1 


3 


4 


1 


1 


5 


1 


2 


6 


1 






E 


rrii 




1 


1 


2 


2 


1 





3 


1 


2 


4 


1 


1 


5 




1 



FP 


mi 




1 


1 


2 


2 


1 





3 


1 


1 



Table 4: Masses und Angular Momenta of the Amplitudes for K — 1 
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